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Abstract 

We present a convergence analysis of the penalty method applied to unilateral contact problems 
in two and three space dimensions. We first consider, under various regularity assumptions on 
the exact solution to the unilateral contact problem, the convergence of the continuous penalty 
solution as the penalty parameter e vanishes. Then, the analysis of the finite element discretized 
penalty method is carried out. Denoting by h the discretization parameter, we show that the 
error terms we consider give the same estimates as in the case of the constrained problem when 
the penalty parameter is such that e = h. 

Keywords: unilateral contact, variational inequality, finite elements, penalty method, a priori 
error estimates. 

AMS Subject Classification: 65N12, 65N30, 35J86, 74M15. 



1. Introduction 

The penalty method is a classical and widespread method for the numerical treatment of con- 
strained problems, in particular the unilateral contact problems arising in mechanics of de- 
formable bodies which involve a nonlinear boundary condition written as an inequality (see, 



e.g., [13, |13, l25| ) . Nevertheless, and to the best of our knowledge, the convergence analysis of 
the method in the simplest case of linear elastostatics with or without finite element discretiza- 
tion has been object of few studies. We may nevertheless quote the earlier, and pioneering works 



of Kikuchi, Kim, Oden and Song [lij2i,[3 (see also Q) and the more recent study dealing 



with the boundary element method 
In the context of elliptic partial differential equations, the penalty method is also classical for 
the treatment of Dirichlet boundary conditions, and has been thoroughly analyzed for instance 
in However, due to the very different nature of Dirichlet boundary conditions and contact 
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conditions, and of the resulting penalty methods, the aforementioned analysis can be hardly 
adapted to contact problems. 

In this paper, we present a convergence analysis of the penalty method for unilateral contact 



which use in particular some recent results from [15|]. We analyze both the continuous and 
discrete problems. We limit the analysis to a conformal discretization with piecewise linear 
finite elements. In particular, we show that the same (sub-optimal, quasi-optimal or optimal) 
convergence rates as for the constrained problem (governed by a variational inequality) can be 
recovered with the choice e = h, where e is the penalty parameter, and h is the mesh size. A 
remarkable fact is that this choice is independent of the regularity of the continuous solution. 
Let us introduce first some useful notations. In what follows, bold letters like u,v, indicate 
vector or tensor valued quantities, while the capital ones (e.g., V, K, ...) represent functional 
sets involving vector fields. As usual, we denote by (H s (-)) d , s G R,d = 1,2,3 the Sobolev 
spaces in one, two or three space dimensions (see [l]). The Sobolev norm of (H s (D)) d (dual 
norm if s < 0) is denoted by || • \\ St r> and we keep the same notation when d= 1, d = 2 or d = 3. 
The letter C stands for a generic constant, independent of the discretization parameters. 

2. Setting 

2.1. The contact problem 

We consider an elastic body Q, in R rf with d = 2 or d = 3. Small strain assumptions are made, 
as well as plane strain when d = 2. The boundary dVl of O is polygonal or polyhedral and we 
suppose that <9Q consists in three nonover lapping parts Yd, Yjy and the contact boundary Yc, 
with meas(r£)) > and meas(rc) > 0. The contact boundary is supposed to be a straight line 
segment when d = 2 or a polygon when d = 3 to simplify. The normal unit outward vector on 
dQ is denoted n. In its initial stage, the body is in contact on Yq with a rigid foundation and 
we suppose that the unknown final contact zone after deformation will be included into Tc- The 
body is clamped on Yd for the sake of simplicity. It is subjected to volume forces f £ (L 2 (£l)) d 
and to surface loads g € (L 2 (YN)) d - 

The unilateral contact problem in linear elasticity consists in finding the displacement field 
u : — y M. d verifying the equations and conditions ([I])-©: 

diver(u) + f = in 0, 

er(u) = Ae(u) in Q, 

v ' y ' (1) 

u = on Fjj, 

er(u)n = g on IV, 

where cr = (<7jj), 1 < i,j < d, stands for the stress tensor field and div denotes the divergence 
operator of tensor valued functions. The notation e(v) = (Vv+Vv )/2 represents the linearized 
strain tensor field and A is the fourth order symmetric elasticity tensor having the usual uniform 
ellipticity and boundedness property. For any displacement field v and for any density of surface 
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forces er(v)n defined on d£l we adopt the following decomposition in normal and tangential 
components: 

v = -u n n + v t and er(v)n = cr n (v)n + er t (v). 
The nonlinear conditions describing unilateral contact on Tc are: 

u n < 0, cr n (u) < 0, cr n (u) u n = 0, (2) 

and the frictionless condition is 

<r t (u) = 0. (3) 

We introduce the Hilbert space: 

V := |v E ( J ff 1 (^)) d : v = on T D } . 

The convex cone of admissible displacements which satisfy the noninterpenetration on the con- 
tact zone Tc is: 

K := {v e V : v n = v • n < on T c } ■ 

Define 

a(u,v) := / <t(u) : e(v) dO, 
Jn 

L(v) := f fvdn+ I g-vdT, 
Jn Jr N 

for any u and v in V. From the previous assumptions, we know that a(-, •) is bilinear, symmetric, 
V-elliptic and continuous on V x V. We know also that L(-) is a continuous linear form on V. 
The weak formulation of Problem ([H)-©, as a variational inequality (see [l2| and also 1^, 14, 



121), is: 

Find u € K such that: 



a(u, v — u) > L(v — u), Vv £ K. 
Stampacchia's Theorem ensures that Problem @ admits a unique solution. 



(4) 



Lem 



2.2. The penalty formulation of the unilateral contact proble 
Let us first introduce the notation [•]+ for the positive part of a scalar quantity a G 



a if a > 0, 
otherwise. 



In the rest of this paper, we will make an extensive use of the following properties: 



a < [o]+, a[a} + = [a] 2 + , Va G E. (5) 
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We introduce also: 

W = {v n]Fc : v € v| , 

and its topological dual space W , endowed with its usual dual norm. Since Tc is a straight line 
segment or a polygon, we have H^ 2 (Tc) C W C H l / 2 (Tc) which implies W C (if*/ 2 (r i c')) / 
where H^ 2 (Tq) is the space of the restrictions on of functions in H 1 / 2 (dtt) vanishing outside 
Tc and H^^jXc) is the space of the restrictions on Tc of traces on <9f2 of functions in 
We refer to [21[ and [l| for a detailed presentation of trace operators and/or trace spaces. 
Let e > be a small parameter. The penalty method for unilateral contact problem ([J]) reads: 

Find u e E V such that: 

1 (6) 
a(u e , v) + -([u £ .n] +1 v n }r c = L(v), Vv 6 V, 

where (v)r stands for the duality product between W' and W. Note that this formulation 
(which is an approximation of the exact contact conditions ([2])) is obtained by setting the 
condition a n (u £ ) = — ^[u £ ^ n ] + instead of the conditions in ([2]) on the contact boundary Tc- 

Remark 2.1. Note that for Problem it holds : u £! n,v n £ W C L 2 {Tc), and so [u £ ^ n ] + € 
L 2 (Tc)- It results that the duality product {v}r c can also be understood simply as the inner 
product in L 2 (Tc)- This remark still holds for the remaining part of the paper in which all the 
duality products (•, -)r c can be changed with inner products. 

We recall that Problem ([6]) is well-posed using an argument proposed by H. Brezis for M type 
and pseudo- monotone operators 0] (see also 20] and [3]): 



Theorem 2.2. For all e > 0, Problem ([6]) admits a unique solution u £ . 

Proof. Using the Riesz representation theorem, we define a (nonlinear) operator B : V — » V 
with the following formula: 

(Bv,w)i ]S! := a(v,w) + -{[v n ] + , w n ) Tc , Vv,w € V, 

where denotes the inner product in (H l (Q)) d . Note that Problem ([6]) is well-posed if 

and only if B is a one-to-one operator. 

Let v,w G V, it follows from the definition of B that: 

(Bv - Bw, v - w)i )n = a(v - w, v - w) + -(K]+ - [w n ] + ,v n - w n )r c . 

Due to the properties ([5]), we observe that, for all a, b G E: 

([a]+ - [6] + )(a - b) = a[a} + + b[b] + - b[a]+ - a[b} + 
>[a] 2 + + [b] 2 + -2[a] + [b] + 
= ([a]+ - [b] + f > 0. 
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This property combined to the V-ellipticity of a(-, •) imply that there exists a > such that: 
(Bv — Bw, v — w)i n > a||v — w||f jjj Vv,wgV. (7) 

Let us also show that the operator B is hemicontinuous, which means that for all v, w G V, the 
real function 

[0, l]3t^ ip{t) := (B(v - tw), w)i, n € M 
is continuous. For s,i G [0, 1], we have: 

|p(t)-¥>(s)| = |(B(v-tw)-B(v-aw),w) 1)n | 

< |a(v - tw, w) - a(v - sw, w)| + -|([v n - tw n ]+, w n ) Tc - (K - w>n)r c | 



|a((s - t)w, w)| + - 



([v n - tw n ]+ - [v n - sw n ] + ) w n dT 



<\t- s\a(w, w) + - / \[v n - tw n } + - [v n - sw n } + \ \w n \ dT. 
Using then the inequality \[v n ]+ — [w n ]+\ < \v n — w n \, we obtain: 

\<p(t) — f(s)\ <\t — s\a(w, w) H — / \t — s| |w ra | \w n \ dT 

£ he 

= \t - s\ fa(w,w) + ^||^n|lo,r c ^ ■ 

It follows that ip is Lipschitz, so continuous. The operator B is then hemicontinuous. Since 
d?]) also holds, we can apply the Corollary 15 (p. 126) of [7| to conclude that B is a one-to-one 
operator from V to V. That concludes the proof of the theorem. □ 

2.3. Finite element setting and discrete penalty problem 

Let ~V h C V be a family of finite dimensional vector spaces (see 0]) indexed by h coming from 
a family T h of triangulations of the domain O (h = max Te7 -ft where hr is the diameter 
of T). The family of triangulations is supposed regular, i.e., there exists a > such that 
VT € T h ,hx/pT < a where pr denotes the radius of the inscribed ball in T. We choose 
standard continuous and piecewise affine functions, i.e.: 

V h = {v h G (C(Tt)) d : v h | T G (Pi (T)) d ,VT G 7~\ v h = on T D } . 
We introduce W h (T c ), the space of normal traces on Tc for discrete functions in V : 

W h (T c ) := {p h G C(F^) : 3v fe G V h , v h ■ n = 

To simplify, we suppose that the end points (or the border for d = 3) of Tc belong to Tat. 
Moreover we assume that the mesh on Tc induced by T h , on which are defined functions of 
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w h {r c ), is quasi-uniform, which implies in particular that it is locally quasi-uniform, in the 
sense of Bramble & al 0] . 

The discrete version of the penalty method (|6|) for Problem dH) reads: 
Find G V h such that: 

a(u h £ ,v h ) + j([u h e>n ] + ,v*} Tc = L(v ft ), Vv fe G V h , (§) 

where e > is still the small penalty parameter. Using exactly the same argument as in Theorem 
12.21 we see that this problem admits one unique solution. 



3. Convergence analysis of the penalty method 

We present in this section the convergence analysis of both continuous and discrete penalty 
methods. We first state the two main theorems. Next, we give some necessary technical lemmas, 
followed by the proof of each theorem. 

3.1. Main results 

3.1.1. Convergence when e — ¥ 

The penalty formulation introduces a consistency error. The following theorem gives a bound 
for this error as a function of the penalty parameter e. In particular, we recover the well-known 
result that the solution u e of Problem © converges to the solution u of Problem (JH) as e — > 0. 
We suppose that u is more than regular so that all the duality pairings become 1? inner 
products. 

Theorem 3.1. Suppose that u, the solution of Problem belongs to (H^ +u (U)) d with v G 
(0,1/2]. Let u £ be the solution of Problem ([6]). We have the a-priori estimates: 



|u - Uell^n < Ce2 +u \\u\\± +u ^, 



a n (u) + -[u £l „] 



< Ce u \\u\ 



0,17 



< Ce 2u \\u 



(9) 
(10) 

(11) 



with C > a constant, independent of e and u. 



3.1.2. Convergence when e — > 0, h — > 

Next, the following theorem provides the convergence rates for the discrete penalty problem, as 
a function of both discretization parameters e and h. 
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Theorem 3.2. Suppose that the solution u of Problem @ belongs to (i/2 +y (Q)) rf with v € 
(0, 1/2] . The solution of the discrete penalty problem ([8]) satisfies the following error estimates 
in two space dimensions: 



|u-u{N|io + £2 



a n (u) + - [u h e 



o,r c 

i 



+ /l^£2 + /^~2£ u 



/i|ln/i|2 +(/ie)2 +e) || u || 2 



ifO < v < i, 



(12) 



< T n (u) + -K )r 



-f,rr 



< c 



hi +! f +v e-f + + /i 2 "-te3 + /i 2 ""^ llu 



/l2|ln^|2£ 2 +/l|ln/l|2 +(/l£)2 -f £ j 1 1 U 1 1 2 



i/0 < v < |, 



(13) 



iwi/i C > a constant, independent of e, h and u. In three space dimensions, the terms /i2 + 2 +l/ 
(resp. h\\n.h\?) in [Wjl and [T3\) have to be replaced with /i2 + 2 (resp. h±). 

One interesting fact is that the choice e = h in the previous theorem, leads to the same error 
estimates as for the finite element approximation of the variational inequality (see [15()- With 
this choice, the error estimates are straightforward and are given next. 

Corollary 3.3. Suppose that the solution u of Problem belongs to {H^ JrU {Q)) d with v € 
(0,1/2]. Suppose also that the penalty parameter is chosen as e = h. The solution of the 
discrete penalty problem ([8]) satisfies the following error estimates in two space dimensions: 



|u-u^||io+/i2 



o- n (u) + -[< n ] 



< C I 



^+!+- 2 ||u|| 3 



o,r c 
ifO < v< i 



a n {u) + -[u e>n \ + 



-vXc 



(14) 



h\ ln/i|2 ||u|| 2 o 



ifu 



2- 



with C > a constant, independent of h and u. In three space dimensions, the terms /i2 + 2 +l/2 
(resp. h\ ln/i|2 ) in |7^p have to be replaced with h^ + 2 (resp. ). 

Remark 3.4. Note that the quasi- optimality is only due to the estimation of the contact term, 
and is also present for standard finite element discretization |7j /. This is not an intrinsic 
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property of the penalty method. Note that optimality could be recovered in both 2D and 3D cases 
(and for any < v < 1/2) ifadditional assumptions on the transition area between contact and 
non contact are made (see fla])- 

Remark 3.5. If u belongs to (H s (n)) 2 with s > 2 (U C R 2 ) then it is easy to show (by using 
the result in 1 24] ) that the error term in ( |i^| ) is bounded by C7i||u|| S) n. 

Remark 3.6. To ensure convergence of the term ||u — u^U^n + ||<r n (u) + i[M^ n ]+|| , 

we need to choose at least e < Ch" v (see fHfy)- As stated previously, to recover optimal 
convergence rate of the terms IW) and IWi) . the best choice is e = h. 

3.2. Preliminary technical lemmas 

Let V h : L 2 (T C ) W h (T c ) denote the L 2 (r c )-projection operator onto W h (T c ). We recall in 
this lemma the stability and interpolation properties of V h : 

Lemma 3.7. Suppose that the mesh associated to W h (Tc) is locally quasi-uniform. For all 
s E [0, 1] and all v € H s (Tc), we have the stability estimate: 

\\V h v\\ s>rc <C\\v\\ s , rc . (15) 
The following interpolation estimate also holds: 

\\v-V h v\\ ,r c <Ch s \\v\\ s ,r c , (16) 
for all v E H s (Tc)- The constant C > is in both cases independent of v and h. 
Proof. The stability estimate (|15p is proven in [6J]. The interpolation estimate comes from e.g. 

8- □ 

We need another lemma which concerns the existence of a discrete bounded lifting from Tq to 

n-. 

Lemma 3.8. Suppose that the mesh on the contact boundary Tc is quasi-uniform. There exists 
U h : W h {T c ) -> V h and C > 0, such that: 



K h {v h )\ Vc -n = v\ \\n h {v h )\\ x> n < C|Klli,r c > ( 17 ) 

for all v h G W h (T c ). 



Proof. The existence of such an operator is proven in [5( (see also [10(] ) . □ 
The proof of each theorem relies strongly on appropriate estimates for the approximate contact 
condition on Tq in dual norm. We give these estimates in the following lemma: 
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Lemma 3.9. Suppose that u, the solution of Problem belongs to (H 2 +u (Q^ d with v € 
(0,1/2]. Let u e be the solution of Problem (|6|). We have the bound: 



<7 n (u) + ~[u E>n ] + 



< C e u 



<7 n (ll) + - [u £ . n } + 



o,rv 



+ e" _ s Hu-UellLn , (18) 



mitt C > a constant, independent of e, u and u £ . 

Lei &e i/ie solution of Problem ^ . We have also the bound: 



o-„(u) + -[«*]+ 



< C \ h v 



O-n(u) + -[«*] 



+ /i^2 ||u-u*||io 



(19) 



o,r c 



wi/i C > a constant, independent of e, h, u and u^. 



Proof. We first prove the bound (|19|) for the discrete solution. The bound (|18p for the continuous 



solution will be derived in a quite similar manner. By definition of the norm 
the projection operator V h : 



- v r c and using 



a n (u) + -[u*]+ 



-v,r G 



sup 



(a n (u) + - £ [u h £ J+,v) Vc 



Mw,r a 



(a n (u) + l[u^ n ] + ,v-V h v)r c , (a n (u) + l[u* n } + ,V h v) Tc 

— SU P Mil SU P j-j — j-j 

ueff"(r c ) IMkr c veH»{r c ) IMkr c 



< 



1, 



<7„(u) + -K ; J + 



sup 



b-^||o,r c t n K(u) + i[«y+,?Mr c 
-r O sup 



<C /i" 



CJ n (u) + J + 



+ sup 

0,r c v£H»(r c ) 



v \\v,T c veH»(Tc) WPhV\\u,T c 

(<r n (u) + ±[u*]+,V h v)r c \ 



\\V h v\ 



v,Tc 



(20) 



In the fourth line, we used the Cauchy-Schwarz inequality and the stability property (|15h . The 
last line is a direct consequence of the interpolation property ()16j) . Now, with help of the 
relationship a(u — u £ ,v h ) = (<r n (u) + n ] + , u£)r c for all v h G V fe , of the discrete lifting (fT71) 
and using the continuity of a(-, •), it results that: 
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sup 



(a n (u) + i[u*]+,V h v)r c (a n (u) + ±[u*]+,K h {V h v)\r c • n)r c 



\V h v\\ u ,r c 



sup 
sup 



WP h v\\v,To 

a(yL-u h £ ,n h {V h v)) 



\\V h v\ 



<C||u-u£||i, n sup 



*ff"(r c ) 



IPS! 



<C||u-u £ ||i,q sup 



Since the discrete trace space W ft (rc) has the quasi-uniform mesh property, we make use of the 
inverse inequality 

\\V h v\\i Tr <Ch v -^\\V h v\\ v , Tc , 



so that we finally get: 



sup < Ch" 2 u -u^ i, n . 

vZH"{T c ) \\P v \\v,Tc 



This, together with (|20|) leads to: 



0n(u) + -[t£j- 



o- n (u) + -[«; 



h 1 

£,nJ 



u — u: 



o,r c 



which is the desired bound (fT9j) . 

For the continuous problem, we introduce V £ , a fictitious finite element space, defined identically 
as ~V h and with the choice of mesh size h = e. We define also a fictitious discrete trace space 
W e (Tc), in the same manner as W h (Fc), and with also the mesh size h = e. We note simply 
V s : L 2 (Tc) — > W £ (Fc) the L 2 (rc)-projection operator on W £ (Tc)- We again obtain an 
analogous estimate as in (j20|) : 



a n (u) + ~[u £ ,n] 



<C\e v 



cr n (u) + - [u e ,n]- 



+ sup 

0,r c v&H v (T c ) 



(a n (u) + ±lu £ , n ] + ,V £ v)r c 

WP £ v\\v,T c 



With help of the relationship 



a(u - u e , v) = (cr n (u) + ~[u £jn ] + , v n )r , v € V, 



(21) 
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using the continuity of a(-, •) and a discrete lifting operator 1Z £ , we come to the conclusion that: 



0"n(u) + -K, n ] 



< C \e v 



+ e 2 u - u e i o 



o,r c 



which is the desired result (fTHI). 



□ 



Remark 3.10. When v = 1/2 i/ie previous result can be improved, by using simply estimate 
HH) to obtain 



<7 n (u) + -[w e , n ] 



-i/2,r c 



< C||u - u e ||i,n. 



5.5. Proof o/ Theorem\3J\ 

We first use the V-ellipticity of a(-, •) and the fact that u (resp. u £ ) is the solution of Problem 
flD (resp. ©) to obtain: 



«ll u - u e|ll,n < a ( u - u e, u ~ u e) 

= a(u, u - lie) - a(u e , u - u e ) 

= (o"n(u) + -[«£,„]+,«„ - -U £ ,n)r c 

= (f7 n (u),U n )r c + (^[u £i „] + ,u n )r c - (cr„(u),u £in )r - {^[Ue,n]+,Ue,n)T c , 
where a > denotes the ellipticity constant. Contact conditions ([2]) on Tc yield: 

(a n (u),u n )r c = 0, (i[n £in ] + ,u n )r c < 0. 
Once again conditions ([2]) associated to properties ([5]) provide: 

-(a n (u),u e>n )r c < -(cr„(u), [n £jn ] + )r c , 

{-[u £ ,n] + ,U £ ,n)r c = (-[V>e,n] + , K,n]+)r C - 

We then get: 



«ll u - u e|ll,fi < ~(<7n(u) + -[« e ,n] + > K,n] + )r c - 



So we continue bounding: 
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«ll u - u elli,n 

< - (e(a n (u) + ~[u e>n ] + ), ^[n £jn ] + + er n (u) - cr n (u))r c 

2 



<-e 



<-e 



1, 



cr n (u) + -[tl e , n ] 



cr n (u) + -[%,„]- 



£7 n (u) + -[l%, n ]- 



+ e(cr n (u) + i[u ein ] + ,cr n (u))r c 



o,r c 
2 

o,r c 
o,r c + 2/3 



CJ n (u) + -[lt e , n ]- 
CJ n (u) + -[M e ,„] + 



e 1 *||an(u)||„^ 



*/,r c 

2 



+ 



-v,rc 



kn(u)L, 



r > 



with 5 G [0, 1], P > 0. Note that since we supposed u G (F§ + ^(0)) 2 , we have cr n (u) G H V {T C )- 
We combine this result with estimation (|18|) : 



«ll u - u elli,n <-£\l-C- 



P 



cr n (u) + -[tl e , n ]- 



o,r c 



2(5+//)— 1 



p 



u - u 



+ 



Pe 



2-28 



kn(u)||J ro , 



which can be transformed into: 



< 



/9e 



C- 

2-2<5 



2(5+i/)-l N 
||<7 n (u >M - 



U 



U £ ||?o+£ 1-C- 



-2(5+i/) -1 





<7n(u) + 



o,r c 



Taking 5 = 1/2 — u, P = 2Cmax(l,a 1 ) and using the estimate ||o*n(u)||i/,r — C|l u l 
proves the bounds @ and (|10p of the theorem. The bound ([lip of the theorem is then a direct 
consequence of this last result and the estimation in norm of the negative exponent Sobolev 
space (JIB]). □ 



3.^. Proof o/ Theorem\KE 

We denote by Z' 1 the Lagrange interpolation operator mapping onto ~V h . We first use the 
V-ellipticity and the continuity of a(-, •), as well as Young's inequality, to obtain: 
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«l|u-<lll,n < a(u-<,u-<) 

= o(u - u£, (u - l h u) + (X h u - u*)) 

< C||u - u^||i,n||u - I h u||i, n + o(u - u^,I ft u - u£) 

< f llu - u^||?, n + ^llu -^ull?,,, + a(u,I h u - u h e ) - a(u h £ ,l h u - u h e ), 

with a > the ellipticity constant. Since u is solution of ([!]) and is solution of (jHJ), we can 
transform the term a(u,X h u — u^) — a(u^,I h u — u^). So we obtain: 



-||u - u^\\l n < — ||u - l h u\\l n + <<r n (u) + -[<„]+, (I h u) n - uJ n ) ro . 

Because of conditions ([3j) and since the lD-Lagrange interpolation with piecewise-linear poly- 
nomials preserves the positivity, we have that (I h u) n < on Tc- This implies: 

<~[<J+,(2*u) n >r o <0. 
Once again condition ([3]) associated to properties ([5]) yield: 

-(a n (u),u* n )r c < -((T n (u),[u^ n ] + )r c , 
(i[<J+,< n )r c = (^[n^ n ]+,[u e h n ] + )r c . 

This results into: 



ol C 1 

f llu - < 2^[|u - 2*u||? in + ( ff „(u), (X h u) n )r c - (a n (u) + -[<„]+, [<J+)r c . (22) 

We bound the last term (a n (u) + -[it^ J + , [u^ n ] + )r c as follows: 



(o-n(u) + -K, n ] +> [«* n ]+)r 

(e(cr n (u) + + <r n (u) - cr n .(u))r c 

2 



<-£ 



<Tn(u) + -[«J n ]- 



^n(u) + -[<„]- 



+ e(<r n (u) + -[«£„]+, 0Vi(u))r o 
o,r c 

2 



o,r c 



o- n (u) + -[uj 



<T„ U 
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The last inequality has been made possible due to the assumption on the regularity of u, which 
implies that cr n (u) £ H u {Tc)- Treating the last term with Young's inequality, and with (3 > 0, 
we obtain finally: 



1 



(a n {u) + -[<„]+, K,J+)r c 

1 2 
1 ' h 



< -e 



c%(u) 



o,r c + ^ 



a n (u) + -[u h £ ] 



+ 7rlkn(u)||^ rc . 



(23) 



We now combine this last inequality (|23p with (|22p and then insert the estimation (|19p into the 
resulting inequality: 



"'In ii^ll 2 

l U U £ II 1,0 



2 

<^||u-T"u|jio + (.7, i (u).(X ,i u)„} l . v 



■-fa. 1 1 2 

ll,: 

1 



<Tn(u) + -[<J 



o,r„ + 2/3 



(7 n (u) + -[<«] + 



-vXc 



+ f lkn(u)|| 2 >rc 



c 2 



£ 1-C 



fe..ll2 

1 

eft 2 " 



< ^-||u - l n u\\( >n + (a n (u), (^u),,)^ 



2/3 



ow(u) + -[«*,„]- 



o,r c 



+ C/ i 2 - 1 ^||u-u £ || 2 if , + ^K(u)|| 2 irc . 



We rearrange the terms: 



C/^" 1 ^- ||u-u e h || 2 n +£ 1-C 



eh 2u 



o-„(u) + - [u: 



h ] 

e.ni 



< ^-\\ u - l h u\\l n + K,(u), (Z fc u) n ) rc + ^K(u)|| 2 rc . 



We choose /3 as follows: 



o,r c 



(24) 



P = C {h 2v - x e 2 a x + h 2 »e) , 

with C > a constant sufficiently large so that the two left terms in (|24|) are positive irrespec- 
tively of the values of e and h. It results that: 



/3K(u)|| 2 jrc = C(h 2 »- 1 e? + h*'e) K(u)|| 2 jrc . (25) 

The estimation of the Lagrange interpolation error in L 2 and H 1 norms on a domain D is 
classical (see e.g., [Ill]): 
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/i- 1 ||u-X /l u|| 0i D + ||u- J ft u||i )D < Ch^WuW^r,, (26) 

for s € (1,2]. The contact term (<r n (u), (I h u) n )r c can be estimated in two space dimensions 
using results from [151 ]: 



^+,+2^|| u ||2 



Wu),(lH)r c <C , 



2 

2,n 



if < v < 5, 



(27) 



In three space dimensions the bound is obtained in a straightforward way using ([26]) for any 
< v < 1/2 : 



Mu),(2 fc u) fl > ro <Crfc 1+, '||u|| 2 |+| , >n . 
In two space dimensions, we combine finally the estimations (|24p — (]2T|) to prove that: 



(28) 



|u-«£lli,n + E 



<7„(U) + -K, n ] + 



o,r c 



l+u+2u 2 , f,2i/-l_2 , 1,21/ 



if < v < i 



2' 



(/i 2 | ln/t| + e 2 + /te) ||u| 



2 

2,fi 



if l/ = i, 



which is the first required estimate. This, together with the estimate (|19p . yields additionally 
the bound in two space dimensions: 



0W(u) + -[uJ n ] 



< C ^ 



((fc"e _ 3 + ^-§)/ t l+f+ i/2 + ^-I e l + ^ + /i 2 ^) ||u||| +I/)n if < v < ±, 
(h2 £ -2 + l)h\ ln^|5 + h2 £ 2 +h + e) [|u|| 2 ,n iff — 1 



2' 



Using /i 2l/ < /i 2 +!y and /i < h\lnh\2 ends the proof of (I13p . The bounds in three space 
dimensions are obtained as before using estimate (|28|) instead of (|27p . □ 
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